Phonon Universal Transmission Fluctuations and Localization in Semiconductor 
Superlattices with a Controlled Degree of Order 
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r^i , We study both analytically and numerically phonon transmission fluctuations and localization 

in partially ordered superlattices with correlations among neighboring layers. In order to generate 
a sequence of layers with a varying degree of order we employ a model proposed by Hendricks and 
Teller as well as partially ordered versions of deterministic aperiodic superlattices. By changing 
' a parameter measuring the correlation among adjacent layers, the Hendricks- Teller superlattice 

exhibits a transition from periodic ordering, with alternating layers, to the phase separated opposite 
limit; including many intermediate arrangements and the completely random case. In the partially 
ordered versions of deterministic superlattices, there is short-range order (among any N consecutive 
^ , layers) and long range disorder, as in the N-state Markov chains. The average and fluctuations in 

the transmission, the backscattering rate, and the localization length in these multilayered systems 
are calculated based on the superlattice structure factors we derive analytically. The standard 
deviation of the transmission versus the average transmission lies on a universal curve irrespective 
of the specific type of disorder of the SL. We illustrate these general results by applying them 
to several GaAs-AlAs superlattices for the proposed experimental observation of phonon universal 
transmission fluctuations. 

' 62.65.+k, 63.50.+X, 68.65.+g, 71.55.Jv 
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I. INTRODUCTION 



, For a long time, electronic devices were made of a sin- 
OS ' gle semiconductor material. This is no longer the case. 
Epitaxy and heterostructures have brought a revolution 
in device technology by placing different semiconductors, 
| with different physical properties (dielectric constants, 
energy gaps, etc.), within distances of a few nanome- 
| ters. For example, different lattice sizes in different ad- 
jacent semiconductors produce strain in the heteroepi- 
taxy, altering its physical properties. Furthermore, re- 
cent developments in the technology for stacking different 
semiconductors, in order to fabricate multilayered thin- 
films, makes possible the realization of various semicon- 
ducting superlattices (SL's) with artificially imposed one- 
dimensional (ID) order in the growth direction. Specifi- 
cally, in addition to the usual periodic stacking of semi- 
conductors, several aperiodic multilayers have been fab- 
ricated, including quasicrystalline, Thue-Morse, and ran- 
dom superlattices. Their physical properties have been 
studied by a variety of experimental probes, including X- 
ray and Raman scattering. For a review on these topics, 
with further references, the reader is referred to Ref. [|. 
We note that the experimental studies of acoustic wave 
propagations (of both phonon and ultrasonic regimes ) 
in some of these aneriodic systems have also been done 
by several groupsB'H. 

It is the purpose of this work to study systematically 
the phonon transport properties of superlattices as a 
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function of their structural order. In particular, we study, 
both analytically and numerically, the transmission fluc- 
tuations and localization properties of phonons in two 
types of partially ordered SL's, which are described in 
more detail in the next two sections. The .first type 
is based on the Hendricks- Teller (HT) modelo for lay- 
ered systems, which has the very convenient feature of 
having a tunable degree of structural correlation among 
neighboring layers. In particular, we consider the grad- 
ual and systematic transition from a periodic arrange- 
ment of alternating layers to the opposite, phase sepa- 
rated, regime and follow the corresponding changes in 
the transport properties induced by the changing struc- 
tural order of the SL. The second type is based-on the so- 
called three- and four-state Markov structuresB, and is il- 
lustrated with two examples, which am.partially-ordered 
versionsti o£the quasicrystalline (QC)Lro and Thue-Morse 
(TM) SL'stI It should be noted that the random version 
of QC SL's as defined by a three-state Markov process 
were also fabricated and the Raman spectra in these sys- 
tems have already been measuredH. 

Our strategy is the following: we derive analytical ex- 
pressions of J s , the average phonon intensity reflected 
from the interface of layers, for several SL's with a vary- 
ing degree of short-range correlations. From / s , we 
analytically derive the localization lengths, transmission 
rate, and transmission fluctuations, all of which coincide 
well with the numerical results we obtain from the al- 
ternative transfer matrix method. The relation between 
the different quantities which characterize phonon trans- 
port is presented. For instance, the Lyapunov exponent, 
which provides the inverse of the phonon localization 
length, is the logarithmic decrement of the transmission 
coefficient averaged over the realizations of disorder. We 
apply the general ideas and results derived here to sev- 
eral particular realizations of GaAs-AlAs SL's which are 
readily accessible experimentally. Our predictions for 
the universal phonon transmission fluctuations can be 
tested using currently existing experimental techniques 
in phonon spectroscopy and phonon imaging which have 
so far been used to verify the existe±ice-Qf phonon filter- 
ing actions of periodic and QC SL'sBclii3. The analogies 
and differences with the universal conductance fluctua- 
tions for transport in disordered systemscJ and speckle 
phenomenally will also be discussed. 

In sections II and III, we describe in detail the two 
novel families of superlattices considered here. In section 
IV, the phonon backscattering rate is studied in terms 
of the structure factors of the SL's, which we derive in 
closed form. Section V is devoted to the average trans- 
mission, transmission fluctuations and the Lyapunov ex- 
ponent. We illustrate in Sec. VI our general results by 
applying them to several proposed GaAs-AlAs superlat- 
tices, for the experimental observation of phonon trans- 
mission fluctuations. Section VII presents a summary of 
our results. 
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II. HENDRICKS-TELLER SUPERLATTICES 



Consider a SL with two kinds of layers, hereafter de- 
noted by A and B, occurring with frequencies /a and 
Jb {Ia + Ib = 1, Ia > Ib)- To introduce a correla- 
tion, consider two adjacent layers and denote by Qaa 
the probability that layer A is followed by layer A, Qab 
the probability that A is followed by B, and so on. The 
first layer of the pair is A (B) with probability Ja (Jb), 
thus 

Qaa + Qab = /a, Qba + Qbb — /s; (1) 
and similarly 

Qaa + Qba = /a, Qab + Qbb = Jb ■ (2) 

From Eqs. (|l|) and ^ we find 

Qaa = /a - 1/4 + q, Qbb = Jb - 1/4 + q, Qab = Qba = 1/4 - q, 

(3) 

where 1/4 — f a < q < 1/4, and q measures the degree 
of correlation among adjacent layers. Note that for Ja = 
/b = 1/2, q = (i.e., no correlation) corresponds to 
the completely disordered case. Let Pab = Qab 1 1 a be 
the probability that the second layer of the pair is B if 
layer A is now introduced as the first layer of the pair. 
This conditional probability can be defined for any pair 
of layers (e.g., Pba)- For convenience, we also introduce 
P(aa\b) describing the conditional probability that the 
B layer is generated after the pair of layers AA, and so 

In the Hendricks-Teller structure!, the probability that 
a layer is present in a certain position depends on the 
neighboring layers as well as the abundance of the layer 
in question. By changing the value of the parameter q 
measuring the correlation among neighboring layers, it 
is possible to conveniently obtain a variety of different 
arrangements ranging from the alternating checkerboard- 
like periodic pattern to the phase segregated case. For 
a positive value of q, the same kind of layers tend to 
stack side by side; thus, as q increases, the number of 
the interfaces between the layers A and B decreases. For 
q = 1/4, the system becomes a phase separated single 
hetero-structure where every layer of material A (B) is 
attached to material A (B), except at the only interface. 
The case q = corresponds to the completely random SL, 
which has been studied in detail in Ref. [l||. For a negative 
value of q, layers A and B tend to stack in an alternating 
fashion, and as q decreases the system becomes closer 
to a periodic SL which is attained for q = —1/4 ( with 
Ja = Ib = 1/2). In summary, a negative q encourages 
alternation among layers while a positive q favors phase 
segregation. 

The study of quasicrystalline diffraction patterns has 
been partly, responsible for a renewal of interest in the 
HT modelli3. Several variations of it have been consid- 
ered. In one of them, the independent random variables 
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are the spacings between the planes (or scatterers). In 
another one, the planes are first periodically spaced and 
then randomly displaced. This difference is not-ohysi- 
cally significant from the point of view of spccklcolij. 

III. N-STATE MARKOV SUPERLATTICES 

Let us now consider a different type of SL with a con- 
trolled degree of randomness. It is modelled after the 
so-called Markov property in thje theory of fluctuations, 
noise, and stochastic processes^. In fact, the subclass of 
Markov systems is by far the most-important stochastic 
process in physics and chemistry3'La Since SL's based 
on this structure are not well known in the multilayer 
community, it is worthwhile to explain the origin and 
motivation for this kind of system, and a few results use- 
ful for calculations in the next few sections. 

The oldest and best known example of a Markov pro- 
cess in physics is Brownian motion. If a series of obser- 
vations of the same Brownian particle gives a sequence 
of locations, fx, r*2, r n , f n +i, ... , each displace- 
ment, <5„+i = r n +\ — r n is affected by chance, and its 
probability distribution only depends on r n , and is inde- 
pendent of the previous history r> t „2, ■ ■ ■ • Thus, 
on the sequence of time intervals imposed by a particular 
experiment, the position and the velocity of the particle 
are Markov processes. This picture forms the basis of the 
theory of Brownian motion. Other examples of Markov 
processes are: the radioactive nuclear decay, the escape 
of gas molecules through a small leak, the destruction of 
cells by radiation, and the emission of light by excited 
atoms. In all of them, on the sequence of time intervals 
imposed by a series of measurements, the state of the sys- 
tem at time t n , only depends on the state at time i n _i, 
and is independent of the states at all previous times, 
t n -2, t n -3, ■ ■ ■ ■ Also, the concept of a Markov process 
is not restricted to one-component processes, but applies 
to m components as well. The three velocity components 
of a Brownian particle and the m chemical components 
of a reacting mixture are two examples. 

In the previous paragraphs, the word 'process' was 
used in its standard physics manner, i.e., usually referring 
to the time evolution of a system. However, sometimes 
the underlying evolving variable is not time, but space. 
An example is given by the following Markov process: 
the loss of cosmic ray electrons in an absorbing material, 
the traversed thickness playing the role usually assigned 
to time. Here, we also consider space (along the growth 
direction) in a similar manner. Also, in this paper as in 
most experiments, we build our structures one layer at a 
time. Starting with layer X, we add the next layer, either 
X or Y, according to the probabilities Qxx and Qxy- 
Thus, in a Markov SL the addition of any new layer only 
depends on the type of layer (or block of layers) most 
recently added, and not on the previous ones. 

We now proceed to describe Markov SL's with short- 
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range correlations in the sequence of constituent layers. 
Specifically, we consider versions of the quasicrystalline 
and Thue-Morse SL's lacking long-range coherence. The 
deterministic quasicrystalline or Fibonacci sequence has 
long-range order manifested by the presence of a dense 
set of Bragg peaks in its structure factored. It is gen- 
erated by iterating the substitution rules A — >■ AB and 
B — > A, so only three possible neighboring pairs of layers 
AA, AB, and BA (three states) are present, and the BB 
pair never appears. In order to preserve this short-range 
ordering, we generate a Markov sequence based on it, by 
using the following straightforward three-state Markov 
chain rules: (i) layer B is generated with probability one 
after the pair of layers AA , i.e., P(aa\b) = lj because 
AAA is a forbidden arrangement in the original system 
with long range order; (ii) layer A is generated with prob- 
ability one after AB, i.e., P(ab\a) = 1; because BB is not 
allowed in the original structure with long-range order; 
and (Hi) P(ba\A) = T and P(ba\b) = T ~ 2 , respectively, 
where r = (y/E+ 1)/2. The last step is the only one that 
introduces randomness in this structure. Therefore, the 
probability of occurrence (i.e., frequency) for the layers 
themselves are Ja = t _1 and Jb = t~ 2 . All these proba- 
bilities also apply to the original deterministic structure 
with long-range order. However, the Markov sequences 
generated according to the above probabilistic rules lack 
long-range coherence. . 

The Thue-Morse (TM) chainErEj is a deterministic se- 
quence that has a degree of order intermediate between 
the quasiperiodic and random cases. In spite of its ape- 
riodicity, the TM Fourier spectrum exhibits very promi- 
nent peaks that would be absent in a random sequence. 
It is the scaling invariance of the TM chain (periodicity 
on a logarithmic scale) which produces long-range cor- 
relations. Many different prescriptions can generate the 
TM sequence, the simplest one is through the substitu- 
tion rules: A — > AB and B — > BA. In this sequence, the 
adjacent pairs AA, AB, BA, and BB (four states) ap- 
pear with equal probability, and blocks AAA and BBB 
are not allowed. 

We can generate a partially disordered structure, 
which preserves the TM short-range arrangement among 
adjacent layers, by following the simple rules: (i) layer B 
(A) is added after layer AA (BB) with probability one, 
i.e., PtAAlB) = 1 and P(bb\A) — lj add layers A and 
B, with equal probabilities, after the pairs of layers AB 
and BA, e.g., P( A b\a) = 1/2 and P(ba\A) = V 2 - Notc 
that /a = /b = 1/2 holds in the partially ordered TM 
sequences. 

IV. BACKSCATTERING RATE OF PHONONS 

In a recent paper, we have shown that the transmission 
rate, localization length and transmission fluctuations of 
phonons in random SL's are derived from the backscat- 
tering rate of phonons due to mass density fluctuations in 
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Incorporating all of the forward scattering contri- 
butions, we can relate the scattering rate to the ensemble 
average of the squared SL structure factor, 7 S) defined by 
I s = (\S N \ 2 )/N, where 

N j 

S n = E ex P(- ? E ™)> ( 4 ) 

is the structure factor of a SL (9a — 0). Here, N is the 
number of A- and £?-blocks (several consecutive identical 
A-layers define an A-block, or block A) in the SL, and 
N — 1 is the number of interfaces between blocks A and 
B. We use the words "layer" and "block" in the following 
way: an elementary or basic layer made of material A, 
with thickness d,A, is called an A-layer or layer A. Also, 
n consecutive yl-layers form an A-block, or block A. The 
same notation applies to B. In Eq. (Q), 9 rn denotes 
twice the phase factor which phonons gain in passing 
through the mth block of a SL consisting of a disordered 
sequence of layers A and B. More explicitly, starting 
from an A block, #2,3-1 = IkADij-i and 9 2 j = 2k B D 2 j, 
( j = 1,2,... ) where kA and k B are the wave numbers 
of phonons in A and B layers and D 2 j-i and Z?2j are the 
thicknesses of the (2j — l)th and 2jth blocks in a random 
SL consisting of A and B layers, respectively. [ Note 
that D 2 j-i = (nA)2j-idA and D 2j = (113)2 jds, where 
(riA)2j-i and (113)23 are ^ ne number of consecutive A (B) 
layers making the 2j — 1th (2jth) block, and d,A and ds 
are the thicknesses of the basic or elementary A and B 
layers. ] The interface between two consecutive identical 
layers (e.g., the interface between the A and A layers in 
a AA block) does not produce any scattering, therefore, 
the only relevant interfaces are between different types of 
layers or the interfaces between A and B blocks. 

Now, we calculate the intensity I s for the partially or- 
dered SL's. Assuming N is an even number ( N = 2ra), 
we can rewrite Eq. (||) as 

n 2j -I 

S N = E [! - ex P(-*#2j)] exp(-i E 0m)- (5) 

j — 1 rn — 1 

Thus, we obtain 

n 

|S,f=^(l-e^) 

3=1 

n j-l 

-E E (l-e- ie «)(l-e- i92m )e- i(e «- 1+e «- 2+ - +l ' 2 "'+ l) + c.c. (6) 

j=2 m=l 

We consider the case where no correlation exists between 
the thicknesses of the adjacent blocks (this is valid for 
the case we are considering) and put 

(exp(-i6»2j_i)) = (exp(-2ik A D 2 j-i)) = e A , 

(exp(-i9 2 j)) = (exp(-2i/c B £ ) 2j)) = e B - (7) 

Now, it is straightforward to derive the expression of I s 
for \e A e B \ < lM The result is 
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(l-e A )(l-e B ) 



1 - £A£B 



(8) 



To proceed further, we calculate €a and e B for any given 
partially-random sequence of A and B layers. For the 
Hendricks- Teller model, we find the following averaged 
phase factors 



Pab e" 



1 - Paa e" 



Pba e 



-ib 



1 - Pbb e 



(9) 



and averaged block lengths 
Pab 



(D 2J -i) = (D A ) = 



(1 - Paa)'- 



d A , 



(D 2j ) = (D B ) = 
(10) 



where a — 2k Ad a and b = 2ksdB- Similarly, for the 
three-state Markov SL with short-range quasicrystalline 
order we obtain 



1 

e A = -e 

T 



1 



(B 



(11) 



and (D A ) = rd A and (D B ) = d B - For the Markov TM 
SL 



(A 



-2ia 



e~ 2 % 
(12) 



and (D A ) = M A /2 and (D B ) = 3d B /2. 

The explicit expression of I s = 1^ for SL's based on 
HT model (with f A =f B = 1/2) is 



2(1 + 4g)(cos0 - cos<5) 2 



(1 + 4g) 2 (cos</> - cos<5) 2 + (1 - 4<j) 2 sin 



(13) 



where <f> = (a + b)/2 and <5 = (a — b)/2. Also we find I s = 
jM-qc and j s = jM-tm for Markov QC and Markov 
TM SL's as 



jM-QC 



(1 — cosa)(l — cos b) 



cos a + t [2 — cos(a + 6)] — r 2 cos(2a + b) ' 



(14) 



(1 - Pbb)' 



d.B 



tM—tm |(1 ~ cosa)(l — cos 6) (cos a + cos 6 + 5/2) 

s 1 + cos 2 f cos 2 I - I [cos(a + b) + cos(2a + 6) + cos(a + 2b) + cos(2a + 2b)} ' 

(15) 

Now, according to our previous workEl, the elastic 
backscattering rate of phonons due to mass-density fluc- 
tuations in random SL's is given in the Born approxima- 
tion as 

r( w ) = g- r 2 i s , (16) 
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where D = (( D A ) + ( D B >)/2 = L/N is the average 
thickness of one block in the system (L being the total 
length of the SL), and R = (Z A - Z B )/(Z A + Z B ) with 
Zj = piCi ( i = AandB, and Cj is the sound velocity) is 
the amplitude reflection coefficient. Here, we note that 
both the substrate and detector are assumed to be made 
of A material and the homogeneous system consisting of 
only A material is taken as the unperturbed system based 
on which the calculation of the phonon transmission rate 
is developed. Thus, Dq/ca gives the average time for 
phonons to propagate through the length of a single block 
in the unperturbed structure. 



V. PHONON TRANSMISSION RATE, 
LYAPUNOV EXPONENT AND TRANSMISSION 
FLUCTUATIONS 



In Ref. 15 we have derived a formula which relates 
the phonon backscattering rate to the transmission rate. 
Introducing a scaling parameter t = (L/ca)T = L/t ( 
I — ca/T is the elastic mean free path of backscattering), 
the average transmission rate (T) is given by 



(T)= / 
Jo 



27rAtanh ir\ 

TT~\ — ex P 

cosh7rA 



dX . (17) 



This formula was originally derived in the study of 
the electrical conductivity in one-dimensional disordered 
metalsEx ■— ■ 
The Lyapunov exponent 7, defined byc3 

7 = - lim (logT)/2L, (18) 

L — >oo 

is an important quantity which provides the phonon lo- 
calization length £ = j^ 1 . The Lyapunov exponent is the 
logarithmic decrement of the transmission coefficient av- 
eraged over the realizations of disorder. It can be proved 
that 7 = T/(2ca) = so 7 is directly related to the 

structure factor or I s t£l. Also, the standard deviation of 
the transmission 

AT = ((T 2 ) - (T) 2 ) 1 / 2 (19) 

is calculated from Eq. ([l4|) by employing the relation!! 

(T 2 ) = - ffi . (20) 

In the next section, we will present comparisons be- 
tween our analytical and numerical results for a variety 
of SL's. It is important to point out that the plots pre- 
sented below are not fingerprints (or speckle patterns) 
of specific configurations of disorder but averages over 
many realizations of disorder. The term speckle pattern 
refers to the complex interference pattern in the trans- 
mitted intensity as a function of frequency (or the out- 
going direction). Each realization of a random medium 
(i.e., each sample of the statistical ensemble) displays its 
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own pattern, or "fingerprint" , which reflects the specific 
arrangement of the inhomogeneities (e.g., impurities) in 
that sample. This phenomenon, called "speckle pat- 
tern" , is familiar in opticsEJ and it refers to the inten- 
sity pattern formed on a screen by light reflected from 
a rough surface. The detailed study, with experimental 
predictions, of the phonon spectroscopy analog of these 
"fingerprints" will be presented elsewhere. 

The expression 'universal transmission fluctuations' 
clearly does not refer to the phonon analog of 'univer- 
sal conductance fluctuations', but to the fact that differ- 
ent realizations of disorder have fluctuations which fall 
on the very same universal curve for the standard devi- 
ation versus average transmission. In fact, we do ob- 
tain, analytically and numerically (for a variety of SL's), 
a universal curve (AT versus (T)) for the transmission 
fluctuations. Also, universal conductance fluctuations 
are not directly related to localization, while our focus 
here is on localization. Finally, it has been pointed out 
that the notation "universal conductance fluctuations" is 
a misnomer because it refers to a sample-dependent, and 
therefore non-universal, fingerprint. Currently, they are 
more appropriately denoted by the term "reproducible 
conductance fluctuations." 



VI. COMPARISON BETWEEN ANALYTICAL 
AND NUMERICAL RESULTS 

A. Hendricks-Teller model 

Figures 1-5 present calculations for the quantities de- 
scribed above, obtained by using two very different ap- 
proaches. In one of them, we use the analytical expres- 
sions presented in this paper. In the other one, we use the 
transfer matrix method for numerical calculations. In the 
latter method, the displacement and stress fields associ- 
ated with the incident and transmitted waves are con- 
nected each other by the product of the transfer matri- 
ces describing the physical properties of each constituent 
layer of the SL. The transmission rate is expressed in 
terms of elements of the product of the transfer matrices, 
by imposing proper boundary conditions on the incoming 
and outgoing waves. Readers interested in a pedagogi- 
cal introduction to transfer matrices and other related 
techniques, are referred to Ref. |23|. 

In order to verify the accuracy of our predictions, it 
is important to compare the results obtained from these 
two quite different approaches. 

Figures 1(a) to 1(c) plot (T) versus frequency for the 
Hendricks- Teller model with Ja = Jb = 1/2 and for 
q = 1/8, 0, and —1/8. For q = 1/8, the same kind of 
layers tend to stack side by side; q = corresponds to the 
completely random SL; and for q = —1/8, layers A and 
B tend to stack alternatively. For q — l/A the system is 
phase segregated with a single hetero-structure and the 
transmission rate becomes a constant independent of the 
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phonon frequency. For q = —1/4 the system is a peri- 
odic SL and sharp dips in transmission occur due to the 
Bragg reflection of phonons. Figures 1(a) to 1(c) prop- 
erly reflect the features characteristic of these SL systems 
with a highly-controlled degree of disorder. To plot these 
figures we have chosen 34-A-thick GaAs and AlAs as the 
A and B layers, respectively. The average transmission 
reveals various structures including sharp enhancements 
and dips. The former, with (T) ~ 1, are the resonances 
which occur for phonons whose wavelengths match the 
thicknesses of A and B layers. More explicitly, the reso- 
nances occur for cos a = 1 and cos b — 1, or equivalcntly 
at v = i>l n = nci/2di (i = A or B, and n is an inte- 
ger). This can be seen from Eq. ( |l3] ) by noting that 
(cos0 — cos S) 2 = (1 — cosa)(l — cos 6). Numerically, 
the resonance frequencies are = 490 x n GHz and 
v&l = 582 x n GHz. 

The dips in (T) are due to constructive interference 
of backscattered phonons. In the HT SL's the min- 

( B} 

ima of (T) are realized at the frequencies v = Vn = 
n/2(dA/cA + (Ib/cb)- These are the Bragg frequencies 

( B) 

( numerically v n = 266 x n GHz) in the periodic SL's 
consisting of an alternating stacking of A and B layers. 
It should noted that (T) is monotonically decreasing as 
t (oc r ~ Ig) increases, and I s takes its maximum value 
for sin0 = or v = (see Eq. (|l3|)). We find that the 
overall agreement between the analytical and numerical 
results is excellent, even though we observe large fluctua- 
tions in (T) . These fluctuations are small only for T close 
to zero and unity, and remain large even if we increase 
the system size, i.e., the phonon transmission is not a 
self-averaging quantity, as described in Refs. p4|p5[ 

The frequency dependence of the Lyapunov exponent, 
which is proportional to I s , is also plotted in Fig. 2 for 
q = 1/8, 0, and — 1/8 . At the resonances, 7 vanishes 
because I s = and the phonons are delocalized. The 
maximum values jmax of 7 are achieved at the Bragg 
frequencies v„ . At these frequencies I s = 4/(1 + Aq) 
because sin <fi — and from Eq. ( ]l0| ) 

Do = (fAd A + fBd B ), (21) 

1 — 4g 

so Imax ar e proportional to (1 — 4g)/(l + 4q). This ex- 
plains the relative magnitudes of ^ max for different values 
of q shown in Fig. 2. We also note that in the present case, 
the localization lengths are longer than 3000 A, which is 
much larger than the typical wavelength of 35A at a 1- 
THz frequency. Even under this condition, phonons in 
an infinite, partially ordered SL are localized except at 
resonance frequencies-due to the coherent interference of 
backscattered wavescll 

Figure 3 exhibits AT versus (T) plotted together for 
different values of q. It is important to emphasize that 
the data lies on a universal curve irrespective of the 
value of q (magnitude of the correlation) and also of the 
specific type of ordering of the SL. There is a certain 
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amount of fluctuation in the data numerically computed 
via transfer matrices (shown as scattered points in the 
figure). This fluctuation, however, decreases when the 
average is taken over an ensemble consisting of a larger 
number of SL's. We have also explicitly shown (see Fig. 3) 
that the standard deviation AT vanishes for (T) = 1 and 
0, monotonically increases in the range (T) < 0.4, and 
monotonically decreases in the range (T) > 0.4. 

B. QC and TM Markov superlattices 

The average transmission rate versus frequency for the 
Markov versions of the quasicrystalline and TM SL's are 
plotted in Fig. 4 together with the transmission rate of 
the original QC and TM SL's with long-range determin- 
istic order. The basic layers A and B assumed here, are 
the same ones used for the HT model. In the regular 
QC SL considered here, the transmission dips occur at 
frequencies v m ^ n = (m + nr)v / (2r 2 d) , where m and n 
are integers, v is the average sound velocity in the SL, 
and d = rdA + ds ■ When m and n are neighboring Fi- 
bonacci numbers, i.e., (m,n) — (F p -i,F p ) where F p+ i = 
F p + F p _! and (F 0l Fi) = (0, 1), v m , n = v v = T P- 2 v/2d 
holds and a major dip is realized. We have indicated in 
Fig. 4(a) the set of integers (m, n) for several major dips. 
In the regular TM SL we study, large dips in transmis- 
sion happen at frequencies v n and v n /z — f-n/3, where 

v n = Vn ■ We have labeled in Fig. 4(b) the indices of the 
frequencies for several major dips. We see that the small 
dips exhibiting, the self-similar structures characteristic 
of the QC SLo are smeared out in the three-state QC 
Markov SL, producing rather broad transmission dips. 
Similar results can be seen for the four-state TM Markov 
case. 

Figures 5-7 show the results for the three-state QC and 
four-state TM Markov SL's corresponding to Figs. 1-3 of 
H-T model. In Fig. 5, the average transmission rate, (T), 
versus frequency, v, is presented and the agreement be- 
tween the analytically derived results and the numerical 
calculations is excellent. The transmission fluctuations 
are large for an intermediate value of (T) as in the case 
of the HT model. The Lyapunov exponent, 7, versus 
frequency is presented in Fig. 6. The explicit expressions 
of I s (cx 7) for the QC Markov and TM Markov SL's are 
given by Eqs. (|lj) and ([l5]). From these equations we see 
that the resonances (7 = 0) in these systems occur at the 
same frequencies vf^ ( satisfying cos a = 1 or cos 6 = 1) 
as in the HT SL's. Unfortunately, however, we could 
not find any simple explicit analytical expression for the 
frequencies at which the maximum values of 7 are at- 
tained. The standard deviation of the phonon transmis- 
sion rate AT, versus (T) is shown in Fig. 7. In this figure 
the continuous line is the theoretical prediction, which 
agrees well with the numerical points obtained by aver- 
aging over 100 realizations of disorder. Here, it should 
be noted that the analytical results for both (T) and AT 
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are functions of only the scaling parameter t = L/£, the 
system size divided by the mean- free path. Thus, AT 
versus (T) docs not depend on the structures of the SL's 
and the analytical curve in Fig. 7 is identical to that in 
Fig.3 for the HT model. 

VII. CONCLUSIONS 

In summary, we have derived analytical expressions of 
I s , the average phonon intensity reflected from the in- 
terface layers, for several SL's with varying degrees of 
short-range correlation. From I s , we derive the local- 
ization lengths, transmission rate, and transmission fluc- 
tuations, all of which coincide well with the numerical 
results obtained from the transfer matrix method. 

In the SL's based on the HT model, the introduction of 
correlations among neighboring layers drastically changes 
the behavior of the phonon transmission. In particular, 
the rate of transmitted phonons decreases significantly 
with decreasing q. Also, the fluctuations in the average 
transmission are very small close to (T) ~ 1 and (T) ~ 0, 
and become much larger for intermediate values of (T). 

The introduction of disorder in the QC and TM SL's 
produces a decrease in the phonon long-range coherence 
which is reflected in the smearing out of small peaks in 
I s , and cquivalently the smearing out of the small dips 
in the transmission rate. In spite of these quantitative 
differences in the fine structure, the overall qualitative 
behavior is still the same as in the ordered case, in the 
sense that it still exhibits pronounced peaks and dips in 
approximately the same locations as in the original deter- 
ministic SL's. Here we note, however, that the structure 
factors in the original, ordered QC and TM SL's have 
very sharp peaks at v = v p and v = v n / 3 , respectively, 
which means that I s grows in proportion to the system 
size or N at these frequencies. In the Markov SL's, which 
only preserve short-range QC or TM ordering, / s remains 
finite on the entire frequency range of phonons, even if 
the system size is increased indefinitely. 

We have obtained AT versus (T) through two different 
approaches. The data lies on a universal curve irrespec- 
tive of the value of q (magnitude of the correlation) and 
of the specific type of ordering of the SL as demonstrated 
numerically for both HT SL's and two kinds of Markov 
SL's. This is because both (T) and AT are determined 
only by the magnitude of the elastic mean-free-path and 
the system size but does not depend explicitly on the 
details of the structure of random SL's. 

One of the main findings of this work is that the loss 
of long-range order does not produce large qualitative 
changes in the overall structure of the phonon transmis- 
sion rate, while it drastically affects its fine structure. At 
first sight, it might seem surprising to see that the very- 
short-range correlations among neighboring layers dom- 
inate the overall frequency dependent transmission rate 
of the traveling phonons. By increasing the degree of or- 
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dering in the SL in a controlled manner, we find that the 
long-range ordering is responsible for the fine structure 
present in the transmission rate. Furthermore, this effect 
has been precisely quantified through the computation of 
Lyapunov exponents and other quantities useful for de- 
scribing the localized character of the phonons. More- 
over, we have applied these ideas to two families of SL's 
where the degree of order can be systematically changed 
in a convenient manner. Finally, we have applied the 
results derived here to several particular realizations of 
GaAs-AlAs SL's which are readily accessible experimen- 
tally. Thus, our predictions for the phonon universal 
transmission fluctuations can be tested using currently 
existing experimental techniques in phonon spectroscopy. 

Finally, we note that a very interesting experiment 
has recently been done by Kono et al, which studies the 
localization properties of the third sound waves by di- 
rectly measuring the transmission spectra in ID random 
latticesEZI. The observed averaged transmissivity reveals 
the frequency dependence very similar to those given in 
Figs. 1(a) to 1(c), i.e., the periodic oscillation character- 
istic of the presence of both resonances (enhancements in 
transmission) and localizations (dips in transmission ) of 
the waves. 
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FIG. 1. Average phonon transmission rate. (T), ver- 
sus frequency, u, for a Hendricks- Teller SL with 200 layers, 
/a = }b = 1/2, for (a) q = 1/8, (b) q = 0, and (c) q = -1/8. 
These and the following figures consider 34-A-thick GaAs (A) 
and AlAs (B) layers as the basic units for constructing the 
SL's. The continuous line is our analytical result while the 
open circles are obtained from a numerical calculation using 
transfer matrices and averaging over 100 realizations of dis- 
order. The fluctuations in the average transmission are very 
small when (T) ~ 1, and become much larger for the in- 
termediate values of (T), i.e. for 0.2 < (T) < 0.6 . Note 
also that the amount of transmitted phonons decreases signif- 
icantly with decreasing q. This can be understood as follows: 
Similar kinds of layers tend to stack together when q = 1/8, 
thus reducing the number of interfaces with an acoustic mis- 
match, which are the scatterers. For q — —1/8, layers A 
and B tend to stack alternatively increasing the number of 
interfaces (i.e., scatterers). 

FIG. 2. Lyapunov exponent 7 (the inverse of the local- 
ization length £) versus phonon frequency v, for HT SL's 
with q — 1/8 (thin solid line), (dashed line), and —1/8 
(bold solid line). At the resonance frequencies vf^, 7 van- 
ishes, and at the Bragg frequencies , 7 takes its maximum 
value. The localization lengths are about two orders of mag- 
nitude larger than the typical wavelength of 35A at a 1-THz 
frequency. Note that the larger the number of interfaces 
(i.e., for q — —1/8), the sharper the peaks in the Lyapunov 
exponent as described in the text. 

FIG. 3. Standard deviation of the phonon transmission 
rate, AT = ((T 2 ) - (T) 2 ) 1/2 , versus (T), for the same HT 
SL as in the previous figures. The continuous line is the the- 
oretical prediction, obtained from Eqs. (^) and (po|). The 
points are computed using the transfer matrix method and 
an ensemble average over N av = 100 realizations of disorder. 
Fluctuations in the points diminish for increasing values of 
N av . 

FIG. 4. Average transmission rate (T), versus frequency 
v for the Markov versions (bold solid lines) of the (a) QC 
and (b) TM SL's. Also plotted by thin solid lines are the 
transmission rates for the original, deterministic (a) QC and 
(b) TM SL's with 55 and 64 layers, respectively. The physical 
parameters assumed here are the same ones used for the HT 
SL. Several major dips in (T) are labeled (see the text). 

FIG. 5. Average transmission rate (T), versus frequency v 
of (a) QC and (b) TM Markov SL's. The continuous lines 
are the analytical calculations and the open circles are the 
numerical results for the (a) three-state QC and (b) four-state 
TM Markov SL's. The results shown are those averaged over 
an ensemble of 100 random SL's. Each random SL consists 
of 200 basic A and B layers. 



15 



FIG. 6. Analytically calculated Lyapunov exponent, 7, ver- 
sus frequency for the QC Markov (solid line) and TM Markov 
(dashed line) SL's. 

FIG. 7. Standard deviation of the phonon transmission 
rate AT versus average transmission (T). The continuous 
line is the theoretical prediction. Open circles and squares 
are the numerical results for the Markov versions of the QC 
and TM SL's obtained from the transfer matrix method by 
averaging 100 realizations of disorder. 
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